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The mechanism that leads to high-temperature superconductivity in cuprates remains an open 
question despite intense study for nearly two decades. Here, we introduce an SU(4) model for cuprate 
f*) ■ systems having many similarities to dynamical symmetries known to play an important role in 

D ' nuclear structure physics and in elementary particle physics. Analytical solutions in three dynamical 

[j^ \ symmetry limits of this model are found: an SO(4) limit associated with antiferromagnetic order; an 

SU(2) limit that may be interpreted as a d-wave pairing condensate; and an SO (5) limit that may be 
interpreted as a doorway state between the antiferromagnetic order and the superconducting order. 
l ' It is demonstrated that with a slightly broken SO (5) but under constraint of the parent SU(4) 

symmetry, the model is capable of describing the rich physics that is crucial in explaining why 
cuprate systems that are antiferromagnetic Mott insulators at half filling become superconductors 
O ' through hole doping. 

U ■ 

^ ' I. INTRODUCTION 

+J , 

Describing collective motion in a strongly correlated many-body system in terms of single particle degrees of freedom 
is sometimes not feasible. Even if one can solve the problem with the help of large-scale numerical calculations, such 
a practice may not always be interesting because physics can be totally buried in the numerous configurations used 
in the calculation. In the cuprate systems exhibiting high-temperature superconducting (SC) properties, there is a 
substantial point of view (see, for example, Refi) that the many-body correlations are so strong that the dynamics 
can no longer be described meaningfully in terms of individual particles. On the other hand, similar questions were 
also asked long ago in nuclear physics. Nuclear physicists have developed a completely different approach to solve 
the many-body problems. The approach is based on the fact that collective motions in such quantum many-body 
systems are often governed by only a few collective degrees of freedom. Once these degrees of freedom are identified 
and properly incorporated into a model, calculations may become feasible and, most importantly, the physics may 
become transparent. One systematic way of identifying the relevant collective degrees of freedom in a many-body 
system is the method of dynamical symmetries. 

A quantum system exhibiting dynamical symmetries usually contains two or more competing collective modes. In 
nuclear physics, it has long been known that the transition from spherical nuclei, which dominate the beginnings 
and endings of nucleonic shells, to deformed nuclei, which often dominate the middle of shells, is controlled by the 
microscopic competition between long-range quadrupole-quadrupole interactions favoring deformed rotational mode 
and short-range monopole pairing interactions favoring spherical vibrational mode. This competition in nuclear 
physics may be expressed algebraically as a competition between a dynamical symmetry that favors pairing and a 
dynamical symmetry that favors multipole (particle-hole) interactions 2,3 . Besides nuclear physics, the dynamical 
i-^j ■ symmetry method has also been applied in molecular physics^ and particle physics^. 

In cuprates, it is now widely agreed that the mechanism responsible for superconductivity is closely related to the 
unusual antiferromagnetic (AF) insulator properties of their normal states. There are also compelling arguments that 
O ■ the pair mechanism leading to high-temperature superconductivity does not correspond to ordinary BCS s-wave pair- 
^ , ing. Although experimental evidence implicates singlet (hole) pairs as the carriers of the supercurrent, the interaction 
leading to the formation of the singlet pairs appears not to be the traditional lattice phonon mechanism underlying the 
BCS theory, but rather seems to originate in strong electron correlations. Furthermore, the pairing gap is anisotropic, 
with nodes in the k x -k y plane strongly suggestive of d-wave hybridization in the 2-particle wavefunctions. 

Such observations argue strongly for a theory based on continuous symmetries of the dynamical system that is 
capable of describing more sophisticated pairing than found in the simple BCS picture (which is described by a single 
complex order parameter), and capable of unifying SC and AF collective modes and the corresponding phases on an 
equivalent footing. Then such fundamentally different physics as SC order and AF order can emerge from the same 
effective Hamiltonian as concentration variables (e.g., doping parameters) are varied. In this manuscript, we report 
on an SU(4) model that we have recently develop edLSiiLiSi aiming at understanding the mechanism that leads to 
high-temperature superconductivity in cuprates. Theoretical models employing algebras and groups can be found in 
the field of condensed matter (for recent examples, see Refs i 11 ' 12 ). However, to our knowledge, the powerful dynamical 
symmetry method of the type that we shall describe below has not been applied systematically in this field. 
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II. DYNAMICAL SYMMETRY METHOD 

A system possesses dynamical symmetries if the Hamiltonian can be expressed as a polynomial in the Casimir 
invariants of a subgroup chain. For approximately the same period of time that the high-Tc compounds have been 
known, techniques based on dynamical symmetries in fermion degrees of freedom have been in development in the 
field of nuclear structure physics^. There it has proven fruitful to ask the following questions: what are the most 
important collective degrees of freedom in the low-lying spectrum of complex nuclei, what are the microscopic many- 
body quantum operators that create and annihilate these modes, and what is the commutator algebra obeyed by this 
set of operators? 

Systematic investigation of these questions has led to strong confirmation of the following set of conjectures about 
the nuclear many-body system: (1) Low- lying collective modes are in approximate one-to-one correspondence with 
dynamical symmetries in the fermion degrees of freedom. (2) A dynamical symmetry corresponding to low-lying 
collective modes is associated with a Lie algebra and its subalgebras that are formed from a set of fermion operators 
closed under commutation. (3) Different dynamical symmetry subgroup chains arising from the same highest sym- 
metry group are associated with fundamentally different phases of the theory. These phases are unified in the highest 
group. (4) The unification implied by the preceding point suggests that the many low-lying collective states formed by 
systematic filling of valence shells in heavy nuclei are in reality different projections in an abstract multidimensional 
space of the same state. Equivalently, the different states are transformed into each other by the generators of the 
symmetry. Thus, the systematics of collective modes and phase transitions as a function of concentration variables 
are specified by the group structure. 

Motivated by the observation that in a nuclear system at low-energy excitation, valence nucleon pairs tend to 
correlate strongly to form s (spin 0) pairs and d (spin 2) pairs, Wu, Feng, and Guidry introduced s and d fermion pairs 
as the basic building blocks of the Fermion Dynamical Symmetry Model 2 . These pair operators, when supplemented 
with the multipole operators, close either an 50(8) or an Sp(6) algebra, depending on the choice of the basis^. If one 
neglects the fermionic degrees of freedom of these pairs and treats them as bosons, one has the early Interacting Boson 
Model 3 of Arima and Iachello. All the dynamical symmetries emerging from these models have found correspondence 
with the known low-lying collective modes in nuclear systems, such as the spherical vibrational mode, deformed 
rotational mode, and 7-unstable collective mode. 

It has been demonstrated that dynamical symmetries of the type described in Ref^ are realized to remarkably high 
accuracy in the spectrum and the wavefunctions of microscopic calculations using the Projected Shell Modeli2*i». 
Since this model is of the shell-model type based on single-particle degrees of freedom, and is known to give extremely 
good agreement with a broad range of experimental data (see, e.g., Ref»^), this provides rather conclusive proof that 
these dynamical symmetries arc strongly realized in the low-lying states of complex nuclei. This raises the issue of 
whether similar symmetries might be found in other complex many-body fermion systems such as those important in 
condensed matter. 

The dynamical symmetry method applied in cuprates corresponds schematically to the following algorithm: 

0. Assume the following conjecture: All strongly collective modes in fermion (or boson) many-body systems can be 
put into correspondence with a closed algebra defining a dynamical symmetry of the sort described below. This is a 
conjecture, but there is so much evidence in support of it from various fields of physics that it is almost a theorem: 
Strongly correlated motion implies a symmetry of the dynamics described by a Lie algebra in the second-quantized 
operators implementing that motion. 

1. Identify, within a suitable "valence space", degrees of freedom that one believes are physically relevant for the 
problem at hand, guided by phenomenology, theory, and general principles. In the present case, that reduces to 
defining a minimal set of operators that might be important to describe superconductivity and antiferromagnetism. 

2. Try to close a commutation algebra (of manageable dimension) with the second-quantized operators creating 
and annihilating the modes chosen in step 1. If necessary, approximate these operators, or add additional ones to the 
set if the algebra does not close naturally. 

3. Use standard Lie algebra theory to identify relevant subalgebra chains that end in algebras for conservation laws 
that one expects to be obeyed for the problem at hand. In the present example, we require all group chains to end in 
U(l) x SU{2), corresponding to an algebra implementing conservation of charge and spin. 

4. Construct dynamical symmetry Hamiltonians (Hamiltonian that are polynomials in the Casimir invariants of a 
group chain) for each chain. Each such group chain thus defines a wavefunction basis labeled by the eigenvalues of 
chain invariants (the Casimirs and the elements of the Cartan subalgebras), and a Hamiltonian that is diagonal in 
that basis (since it is constructed explicitly from invariants). Thus, the Schoedinger equation is solved analytically 
for each chain, by construction. 

5. Calculate the physical implications of each of these dynamical symmetries by considering the wavefunctions, 
spectra, and transitional matrix elements of physical relevance. This is tractable, because the eigenvalues and eigen- 
vectors were obtained in step 4. Consistency of the symmetry requires that transition operators be related to group 
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generators; otherwise transitions would mix irreducible multiplets and break the symmetry. 

6. If step 5 suggests that one is on the right track (meaning that a wise choice was made in step 1), one can 
write the most general Hamiltonian for the system in the model space, which is just a linear combination of all the 
Hamiltonians for the symmetry group chains. Since the Casimir operators of different group chains do not generally 
commute with each other, a Casimir invariant for one group chain may be a symmetry-breaking term for another 
group chain. Thus the competition between different dynamical symmetries and the corresponding phase transitions 
can be studied. 

7. The symmetry-limit solutions may be used as a starting point for more ambitious calculations that incorporate 
symmetry breaking. Although no longer generally analytical, such more realistic approximations may be solved 
by perturbation theory around the symmetry solutions (which are generally non-perturbative, so this is perturbation 
theory around a non-perturbative minimum), by numerical diagonalization of symmetry breaking terms, or by coherent 
state method as shown later. 

Of course, in practical calculations only a few carefully selected degrees of freedom can be included and the effect 
of the excluded space must be incorporated by renormalized interactions in the truncated space. It follows that the 
validity of such an approach hinges on a wise choice of the collective degrees of freedom and sufficient phenomcnological 
or theoretical information to specify the corresponding effective interactions of the truncated space. 

III. THE SU(4) MODEL 

We now introduce a mathematical formalism that provides a systematic implementation of the dynamical symmetry 
procedure for cuprate systems. 

A. Choice of operators and the algebra 

The basic assumption of the SU(4) model is that the configuration space relevant for competing superconductivity 
and antiferromagnetism is built from coherent pairs formed from two electrons (or holes) centered on adjacent lattice 
sites^i In cuprates the coherent pairs are believed to exhibit d-wave orbital symmetry^, suggesting coexistence of 
two kinds of coherent pairs in a minimal model: the spin-singlet (p) and the spin-triplet (q) pairs. We adopt the 
d-w&ve pairs with structure defined in Refs>ii*i£ as the basic dynamical building blocks of the SU(4) model 

Pu =^f( fc )4T C -fel P12 
fc 

Qij = ^29( k ) c t+Q,i c -k,j m 

k 

where c\ . creates an electron of momentum k and spin projection i,j = l or 2 (=j or j), g(k) = (cosk x — sink y ) is 
the d-w&ve form factor, and Q = (w, w, 7r) is an AF ordering vector. These pair operators, when supplemented with 
operators of particle- hole type Qij and SV, , 

Qij = c \+QA c k ,j Sij — / I c k.i C k,j ~ 2^%> (^) 

k k 

constitute a 16-element operator set that is closed 7 under a U(4) D U(l) x SU(4) algebra if the condition 

g{k) k, sgn (cos k x — cos k y ) 

is imposed. In Eq. J5J, is the maximum number of doped electrons that can form coherent pairs, assuming the 
normal state (at half filling) to be the vacuum. The U(l) factor, Q + = Qn + Q22 = J2k( c t+Qi c k1 + c k+Ql c H)> m 
U(l) x SU(4) is associated with charge-density waves, which is independent of the SU(4) algebra. For no broken 
pairs the charge-density waves are excluded in the symmetry limit and in the following discussion we shall restrict 
attention to the SU(4) group. 

The SU(4) group has three dynamical symmetry group chains^: 



(Pis 



(4) f . 



(1) 



D SO(4) x U(l) D SU(2) S x U(l) 
SU(4) D SO(5) D SU(2) S x U(l) 

D SU(2) P x SU(2) S D SU(2) S x U(l) 



(3) 
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each of which ends in the subgroup SU(2) S x U(l) representing total spin and charge conservation. 
The 15 SU(4) generators are related to more physical operators through the linear combinations: 

= p\ 2 D = pi2 d — wave singlet pairing (2 operators) 

7?^ = ( i 911 2 922 , 9l1 ^ 922 , — i 912 ^ 921 j 7? = (77^ ) ^ d — wave triplet pairing (6 operators) 
q = ^ Qi2+Q2i ^ _ i Qi2-Q 2 i ^ Q11-Q22 j staggered magnetization (3 operators) ( 4 ) 

S = ( Sl2 + Sal , -i Sl2 ~ S21 , Sll ~ S22 ) spin (3 operators) 

n = J]fc,i c l,i c M = + ^ 22 + ^ particle number (1 operator) 

It will also sometimes prove useful to replace the number operator n with the charge operator M, defined through 

m = \{s 11 + s 22 ) = ±{n-n). 

B. The collective subspace 

The group SU(4) has a quadratic Casimir operator 

C au{4) = D*D + ^-tt+Q-Q + S-S + M(M - 4). (5) 

The group is rank-3 and the irreducible representations (irreps) may be labeled by 3 weight-space quantum numbers, 
(cii f2j 03 )• We assume for the simplest implementation of the model a collective c?-wave pair subspace spanned by 
the following vectors: 

|S) = \n x n y n z n s ) = {^{^{^{D^). (6) 
This collective subspace is associated with irreps of the form 

(01,02,03) = (§,0,0), 

where 2f2 is the dimension of the single-particle space in the SU(4) model, and thus Q is the maximum number of pairs 
that can be accommodated in the model space. The corresponding expectation value of the SU(4) Casimir evaluated 
in these irreps is a constant, 

(C S tt(4)) = § (f + 4 ). 

The dimensionality of the full space is 2 . However, the dimensionality of our collective subspace is much smaller, 
scaling approximately as 17 4 : 

Dim(§,0,0) = ^(§ + l)(§ + 2) 2 (§ + 3). 

Thus for small lattices it is possible to enumerate all states of the collective subspace in a simple model where 
observables can be calculated analytically. 



C. SU(4) model Hamiltonian 

The most general 2-body Hamiltonian within the d-wave pair space consists of a linear combination of (quadratic) 
Casimir operators C g for all subgroups g 

g 

where H g are parameters and the Casimir operators C g are2& 



CsO(5) 


— TT ' ■ 7? + S ■ S - 


hM(M-3) 


C*SO(4) 


= Q-Q + S-S 




CsU{2) v 


= D ] D + M{M 


-1) 


CsU{2) 3 


= S-S 




Cu(i) 


= M and M 2 . 
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(<l,S) r\~,S,M) (D]D,M) 

FIG. 1: Dynamical symmetries associated with the (7(4) symmetry. The generators are listed for each subgroup. 



For fixed electron number the terms in M and M 2 in Eq. Q are constant. The term Ho is a quadratic function of 
particle number and may be parameterized as 

Hq = en + iv7i(rt — 1), 

where e and v are the effective single-electron energy and the average two-body interaction in zero-order approximation, 
respectively. Thus the Hamiltonian can be written as 

H = H + V, V = -G D^D -G^ -it - X Q- Q + kS ■ S (8) 

where Go, Gi, \ an d « are the interaction strengths of d-wave singlet pairing, d-wave triplet pairing, staggered 
magnetization, and spin-spin interactions, respectively. Since (C su u\) is a constant, by using Eq. (J3J) we can eliminate 
the 7?^ ■ if term and after renormalizing the interaction strengths the SU(4) Hamiltonian can be expressed as 

H = H' -G{(l-p)D^D+pQ-Q} + K eS S-S, H' = e off n + ±v cff n(n - 1) (9) 

with (1 —p)G = G^ ff , pG = Xeffj and k c s as the effective interaction strengths, and where < p < 1 for the parameter 
p. Since in this paper we primarily address the ground state properties where S = 0, the last term in Eq. Q will 
generally not enter into the later discussion. 



D. Physical interpretation for the dynamical symmetries 

The three subgroup chains of the SU(4) symmetry define three dynamical symmetries with clear physical meanings. 
The symmetry limits SU(2), SO(4), and SO(5), correspond to the choices p = 0, 1, and 1/2, respectively, in Eq. 0. 
The Hamiltonian, eigenfunctions, energy spectrum and the corresponding quantum numbers of these symmetry limits 
can be explicitly given (see Refi). As we shall now discuss, each symmetry limit represents a different possible low- 
energy phase of the SU(4) system. Fig. 1 illustrates these dynamical symmetries and the relation to the corresponding 
phases. 

The pairing gap A (measure of pairing order) and the staggered magnetization (measure of AF order) Q, 

A = G° ff ( J Dt jD )i/2 j Q=( ^4 ) i/2 ) (10) 

may be used to characterize the states in these symmetry limits. We introduce also a relative doping fraction x that 
is related to particle number and lattice degeneracy through 

Since Q — n is the hole number when Q < n, positive x represents the case of hole-doping, with x = corresponding 
to half filling (no doping) and x = 1 to maximal doping. 



1. The SOU) limit 

The dynamical symmetry chain 

517(4) D 50(4) x [7(f) d SU(2) s x {7(1), 

which we shall term the SO(4) limit, corresponds to long-range AF order. This is the symmetry limit of Eq. @ when 
p = 1. The SO(4) subgroup is locally isomorphic to SU(2)p x SU(2)g, where the product group is generated by the 
linear combinations 



F=|(Q + 5), G = I(Q-5) 
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of the original SO (4) generators Q and S. We may interpret the new generators F and G physically by noting that 
if we transform Qij and SV,- defined in Eq. @ to the coordinate lattice space, 

Qij = ^^(~) C ri C rj = c ti c rj ~ c ti c rj 

r r— even r— odd 

r r— even r— odd 

implying that F is the generator of total spin on even sites and G is the generator of total spin on odd sites. Thus, we 
may interpret the SO (4) group as being generated by two independent spin operators: one that is the total spin on 
all sites and one that is the difference in spins on even and odd sites of the spatial lattice. This clearly is an algebraic 
version of the physical picture associated with AF long-range order. 
Furthermore, the ground state has maximal staggered magnetization 

Q = ±0(1 -x) = \n 

and a large energy gap (associated with the correlation Q-Q) 

AE = 2 X off (1 ~ x)Q 

that inhibits electronic excitation and favors magnetic insulator properties at half filling. In addition, the pairing gap 

A = ±G° cS ^x(l-x) 

is small near half filling (x — 0). We thus conclude that these SO(4) states are identified naturally with an AF 
insulating phase of the system. 

2. The SU(2) limit 

The dynamical symmetry chain 

SU{A) D SU(2) P x SU(2) S D SU(2) S x U(l), 

which we shall term the SU(2) limit, corresponds to SC order and is the p = symmetry limit of Eq. ©. In the 
ground state, there exists a large pairing gap 

AE = G ( ^n, 

the pairing correlation is the largest among the three symmetry limits, and the staggered magnetization vanishes in 
the ground state: 

A = \G° cS ^l - x 2 , Q = 0. (12) 
Thus we propose that this state is a pair condensate associated with a c?-wave SC phase of the cuprates. 

3. The SO (5) limit 

The dynamical symmetry chain 

SU(4) D SO(5) D SU(2) S x [7(1), 

which we shall term the SO(5) limit, appears when p = 1/2 in Eq. @. The SO(5) dynamical symmetry has very 
unusual behavior. Although the expectation values of A and Q for ground state in this symmetry limit are the same 
as that of Eq. l|12fl for the SU(2) case, there exists a huge number of states with different number of it pairs that 
can mix easily with the ground state when x is small. In particular, at half filling (x = 0) the ground state is highly 
degenerate with respect to configurations of different number of 7r pairs. The 7r pairs must be responsible for the 
antifcrromagnctism in this phase, since within the model space only 7r pairs carry spin. Thus the ground state in 
this symmetry limit has large-amplitude fluctuation in the AF order (and SC order). This indicates that the SO (5) 
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symmetry limit is associated with phases in which the system is extremely susceptible to fluctuations between AF 
and SC order. 

S. C. Zhang proposed to unify AF and SC states by assembling his order parameters into a five-dimensional 
vector and constructing an SO (5) group that rotates AF order into <i-wave SC order—. The Zhang SO (5) group is 
a subgroup of our SU(4) group, implying that the two models are related to each other. The essential difference 
is that we implement the full quantum dynamics (commutator algebra) of these operators exactly, while in Ref^i, 
the dynamics is implemented in an approximate manner: a subset of ten of the operators acts as a rotation on the 
remaining five operators {D\D,Q}, which are treated phenomenologically as five independent components of an 
order-parameter vector (superspin ft). Thus only 10 of the 15 generators of our SU(4) are treated dynamically in 
the Zhang SO(5) model. If the full quantum dynamics (full commutator algebra) of the 15 operators is taken into 
account, the symmetry for cuprates must be SU(4), not SO(5). 



E. Closure of SU(4) and no-double-occupancy constraint 

Suppression of double occupancy on sites in the copper oxide planes is critical^ in explaining why cuprate systems 
are antiferromagnetic Mott insulators at half filling and become superconductors through hole doping. To realize Mott 
insulator phases at half-filling it is normal to impose a no-double-occupancy rule - the constraint that each lattice 
site cannot have more than one valence electron - by Gutzwiller projection—. We now show an intimate relationship 
between the SU(4) dynamical symmetry and realization of no-double-occupancy constraint in cuprates. 

The relationship can be clearly seen if we express the momentum-space operator sets Q and (|5J , with no approxi- 
mation to the d-wave formfactor g(k), in the coordinate space as 

Pu = { c iAi - °ii c h) 

reA 

Qij ^ ^ ] ( C r,8 C r,j C r,j C r 



(4,^^ - c f,i4,i) ( 13 ) 
Qij = ^ ( C r,i Cr i3 Cf <3 C f : 



E 



Pl2 



reA 



In Eq. (|13|l . lattice site is separated into r = even and r = odd sites. The summation is over even (or odd) sites of 
the lattice, denoted by A, where the ± sign is + (— ) if A is chosen to be even (odd) sites. The quantity Qij is defined 
by Qij = Qij + hSijCl, where cj i (c Tt i) creates (annihilates) an electron of spin i located at r, while c\ i (cf t i) creates 
(annihilates) an electron of spin i at its four neighboring sites, r ± a and rib with equal probabilities (a and b are 
the lattice constants along the x and y directions, respectively, on the copper-oxide plane) 



C r,i 2 ( C r+a,i + C r-a,i C r+b,i C r-b,i 



Explicit commutation operations show that only if the no-double-occupancy constraint is imposed, so that the anti- 
commutator relation 



{c f , 4 ,4,i'} 



is valid, does the operator set (fC^l close under the U(l) x SU(4) Lie algebra^. Thus, the coordinate-space commutation 
algebra of the operators demonstrates that SU(4) symmetry necessarily implies a no-double-occupancy constraint 
in the copper oxide conducting plane. This suggests a fundamental relationship between SU(4) symmetry and Mott- 
insulator normal states at half filling for cuprate superconductors. 

One immediate conclusion^ is that the no-double-occupancy constraint sets an upper limit for the number of doped 
holes if SU(4) is to be preserved: f2 max — jQ e , where f2 e is the total number of lattice sites. Thus, the maximum 
doping fraction consistent with SU(4) symmetry is 

_ O _ 1 

f = fL ~ 4" 
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Beyond this doping fraction the no-double-occupancy condition cannot be ensured and exact SU(4) symmetry cannot 
be preserved. The empirical maximum doping fraction 0.23 ~ 0.27— for cuprate superconductivity may then be taken 
as indirect evidence for a strongly-realized SU(4) symmetry underlying the superconductivity in cuprates. 



IV. COHERENT STATE ANALYSIS 



The soft nature of the SO (5) phase is seen most clearly if we introduce the SU(4) coherent states2I. The result 
of such a coherent state analysis is a set of energy surfaces that represent an approximation to the original theory, 
in which order parameters appear as independent variables. In the general case, these energy surfaces can exhibit 
minima and these minima may appear at non-zero values of the order parameters, implying spontaneous symmetry 
breaking. 



A. The generalized coherent-state method 



The coherent state method is a well-developed theoretical approach to relating a many-body algebraic theory with 
no broken symmetry to an approximation of that theory that exhibits spontaneously broken symmetry. This method 
may be viewed as the extension of Glauber coherent state theory^ (which is built on an SU(2) Lie algebra) to a 
more complex system having an arbitrary Lie algebra structure. It has also been shown to be equivalent to the most 
general Hatree-Fock-Bogoliubov variational method under symmetry constraints, and has been applied extensively in 
various areas of physics and mathematical physics. 

The coherent state \ip) associated with the SU(4) symmetry can be written as 



with the operator T defined by 



T = cxp(?7oopI 2 + V1W12 - n - c -)- 



(14) 



(15) 



In Eq. (2J, |0*) is the physical vacuum (the ground state of the system), the real parameters 7700 and 7710 are symmetry- 
constrained variational parameters, and h. c. means the Hermitian conjugate. Since the variational parameters weight 
the elementary excitation operators p\ 2 and g{ 2 , they represent collective state parameters for a D-tt pair subspace 
truncated under the SU(4) symmetry. 

It is convenient to use a 4-dimensional matrix representation that was introduced in Ref&S. The corresponding 
matrix elements are defined in the 4-dimensional single-particle basis 

{4 T |0*>, 4J0*), c fT |0*), c Fi |0*)}. 

The unitary transformation operator T of Eq. (|14f) may be written in this matrix representation as 



where we define 



X 





-V- 



v+ 




T 



Yi 



Yi X 

-Xt Y 2 



u+ 
u- 



(16) 



u_ 
u+ 



and where unitarity requires 
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The m's and u's, which are related to 7700 and 7710 in Eq. (|15f) . are variational parameters in the matrix representation 
that will be interpreted later as quasiparticle unoccupation and occupation amplitudes, respectively. 

The matrix (|16f) implements a quasiparticle transformation within the D-tt pair space that is further constrained 
to preserve the SU(4) symmetry. The physical vacuum state |0*) is transformed to a quasiparticle vacuum state \tp) 
and the basic fermion operators, 



-0.8 -0.4 O.O 0.4 O-S 

P 

FIG. 2: Coherent state energy surfaces. The energy units are G^Q 2 /4 for figures (a) and (b), and Xcff^ 2 /4 for (c). Ho is 
taken as the energy zero point. Numbers on curves are the effective lattice occupation fractions, with n/Q = 1 corresponding 
to half filling and < n/Q < 1 to finite hole doping. SO(5) symmetry corresponds to p = 1/2 and the allowed range of (3 is 
[— | \J n/Q, i U n/Q] , which depends on n. 



are converted to quasifermion operators, 



{al^al^, a* a* i} with a ri \^)=0, 



through the transformation 



T 



( CrT \ 

V4,y 



|0*> 



(17) 



With above definitions, it is straightforward to write down the corresponding expectation values of one-body and 
two-body terms in Eq. © in the coherent-state representation: 



<Z?t) = (£>) = iO(n + u + +u_«_), (UtD) = (£>t)2 

(4) - (7r,) = -^n(u+t; + -«_«_), (7rt.Tr) = (4)2 

(C> = W+ -v 2 -), (Q-Q) - (S 2 ) 2 

(h) = h(v 2 + +v 2 _), (S-S) = 

M = (Ky) = (S) = (Q X ) = (Qy)=0 



(18) 



where terms that are of order I/O, smaller than the leading terms have been ignored. 



B. Energy surfaces at symmetry limits 

With all the expectation values expressed in the coherent state representation, one can calculate the total energy. 
In Fig. 2, we show the coherent-state energy surfaces as a function of the order parameter (3 for different electron 
occupation fractions n/fi with p = 0, \, and 1. The order parameter j3 is related to the order parameter Q = (Q z ) 

and the electron number n through (Q z ) = 217/?-\/n/2f2 — (3 2 . For p = [SU(2) limit; see Fig. 2(a)], the minimum 
energy occurs at (3 — for all values of n, indicating SC order. For p = 1 [SO(4) limit; see Fig. 2(c)], the opposite 
situation occurs: (3 — is an unstable point and energy minima lie at (3 — ±\\J n/Q, indicating the presence of AF 
order. 

From Fig. 2(b), the SO(5) dynamical symmetry (with p = i) is seen to have extremely interesting behavior: the 
minimum energy occurs at f3 = for all values of n, as in the SU(2) case, but there are large-amplitude fluctuations 
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FIG. 3: Coherent state energy surfaces for slightly perturbed SO(5). The dotted line indicates the location of the ground state 
as n varies and the corresponding f3 value where the energy curves have a minimum. 



in (3. In particular, when n is near O (half filling), the system has an energy surface almost flat for broad ranges of 
j3. This suggests a phase very soft against fluctuations in the order parameters. As n/fi decreases, the fluctuations 
become smaller and the energy surface tends more and more to the SU(2) SC limit. 

A symmetry limit having such a nature is termed a transitional or critical dynamical symmetr^^. A critical 
dynamical symmetry is a dynamical symmetry having eigenstates that vary smoothly with a parameter (usually 
particle-number related) such that the eigenstates approximate one phase of the theory on one end of the parameter 
range and a different phase of the theory at the other end of the parameter range, with eigenstates in between 
exhibiting large softness against fluctuations in the order parameters describing the two phases. Thus, we believe 
that the rich physics in cuprates, in particular in the underdoped regime, must lie near the SO(5) limit, subject to a 
constraint by the parent SU(4) symmetry. 



C. SO(5) symmetry breaking 

Under exact SO(5) symmetry (p = |), the AF and SC states are degenerate at half filling. There is no barrier 
between AF and SC states, and one can fluctuate into the other at zero cost in energy (see the n/tl = 1 curve of Fig. 
2(b)). This situation is inconsistent with Mott insulating behavior at half-filling. The Zhang SO (5) model has been 
challenged because under exact symmetry it does not fully respect the phenomenological requirements of "Mottness" . 
As Zhang— has recognized, for antiferromagnetic insulator properties to exist at half filling, it is necessary to break 
SO(5) symmetry. Such breaking of the SO(5) subgroup symmetry is implicit in the SU(4) model, occurring naturally 
in the SU(4) model if p > 1/2 in the Hamiltonian JjJ. Furthermore, the SU(4) symmetry leads to the following 
constraint 

(£ su4 ) = {D^D + Q ■ Q + t? • tt) = ±fi 2 (l - x 2 ), (19) 

which ensures a doping dependence in the solutions. Thus, the coherent state analysis indicates that the phenomeno- 
logically required SO (5) symmetry breaking and the doping dependence in the solutions occur naturally in the SU(4) 
model. 

To see how in the SU(4) model a broken SO (5) symmetry can interpolate between AF and SC states as particle 
number varies, let us perturb slightly away from the SO(5) limit of p = |. In Fig. 3, SU(4) coherent state results for 
p = 0.52 are shown. One sees from Fig. 3 that if n is near half filling (n/f2 = 1), /3 m ; n = ±0.5; this corresponds to AF 
order, since the staggered magnetization reaches its maximum, Q = 0/2, and there are no pairing. With the onset 
of hole doping, n/fl decreases (x increases). The AF correlation Q quickly diminishes and the pairing correlations 
increase. 

The SO (5) symmetry breaking in the Hamiltonian (p > h) is crucial. Only when SO (5) is broken does the energy 
surface interpolate between AF and SC order as doping is varied (compare the surfaces for p = 1/2 and p — 0.52 in 
Figs. 2(b) and 3). We thus conclude that high temperature superconductivity may be described by a Hamiltonian 
that conserves SU(4) but breaks (explicitly) SO (5) symmetry in a direction favoring AF order over SC order. 
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V. GAP DIAGRAM AT T=0 

We have shown that with the coherent state method, one is able to calculate analytically all the physical quantities 
at and away from the symmetry limits. This provides a simple framework to study in more detail the observables. As 
our first application, we use the preceding results to study energy gaps in the cuprate systems at zero temperature. 

A. The gap equations 

We start from the symmetry-constrained variational Hamiltonian 

H' = H-h\ + G£ sui , (20) 

where H is the Hamiltonian JHJl, A and G are two Lagrange multipliers determined by the constraints of number 
conservation, (h) — n, and conservation of the SU(4) invariant of Eq. I|19|l . 
Introducing the energy gaps 



A d = G' ^(DW) A^G^Tft-T?) A q =x'^/{Q-Q), (21) 

where we define 

Gq = Gq — G G^ = G\ — G x = x — G •> 



one obtains 



G'q G' x x 

Variation of (H') with respect to u± or v± (that is, solving S(H') — 0) yields 

2u±v±{e± - A) - A±(u% - 4) = 0, 



/ A 2 A 2 A 2 \ 
H') = (e-X)n-[ ^ + ^L + -1 . (22) 



which is satisfied by 



with 



e± = V(e±- A ) 2 + A ± A ± = |A rf ±A OT | e±=e T A 9 . (24) 
Inserting Eq. lf2ll|l into Eqs. l(T%)l and (fT^|l . and employing the gap definitions (|2~T]) . one obtains the gap equations 

A d = ^( W+ A+ + W _A_) (25) 

A^ = ^( w+ A+-w_A_) (26) 
4A 

— q - = w + (A q + X')+w_(A q -X') (27) 
X 12 

-2x = w+{A q + A') - w_(A 9 - A') (28) 

l-^ o2 _ fA d \ 2 t fA^ 2 t (A q ^ 2 
4 

where 



"-las + Ur (29) 



u;± = — , A' = A — e 
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By solving the above algebraic equations, all the gaps and the two Lagrange multipliers (the chemical potential A 
and the SU(4) correlation strength G) can be obtained. The three gaps A d , A„ and A q represent, respectively, the 
energy scales of the singlet pairing, triplet pairing, and AF correlation, and we have introduced an analogous SU(4) 
gap Asu4 to represent the energy scale of the SU(4) invariant through 

Asr/4 = G^(£ su4 ) = iGfiVl-z 2 . 
Thus the ground state energy is determined by the four energy gaps (or energy scales) 

E = {H') + nX - \G{1 - x 2 )Q 2 



= ne — 




where (H 1 ) is given by Eq. Ij22(l . Once the gaps and the chemical potential A' are known, the quasiparticle energies 
e±, and the amplitudes u± and v±, can all be determined through Eqs. H23fl - I|24[) . permitting other ground state 
properties to be calculated. 

These results are in many respects analogous to the BCS theory with v± the probability of single particle level e r ± 
being occupied, A± the energy gap, and e± the quasiparticle energies. The essential difference is that conventional 
pairing theories deal with one pairing gap and one kind of quasiparticles; here we have two pairing energy gaps and 
two kinds of interacting quasiparticles, implying a larger variety of possible behavior. The quantities e± are two kinds 
of quasiparticle excitation energies, corresponding to two sets of nondegenerate single particle energy spectra {s±} 
separated by 2A q . Each level can be occupied by only one electron of either spin up or spin down. The corresponding 
pairing gaps are A±, the probabilities of single-particle levels being unoccupied are u±, and the probabilities of 
single-particle levels being occupied are v±. 



B. Solution of gap equations at T=0 



There are three parameters, %, Go, and Gi, in the coupled algebraic equations H25|) - (|29[) . corresponding to the three 
elementary interactions in the SU(4) model: the AF correlation (x), the spin-singlet pairing (Go), and the spin-triplet 
pairing (Gi). Physical solutions of the gap equations depend on the choices for these parameters. Experimental 
evidence suggests that these three interactions in the cuprate system are all attractive, and as demonstrated in RefpiS, 
X should be the strongest and Gi the weakest. Detailed analysis indicates^ that SU(4) symmetric solutions exist 
only when G — G\ (that is, when G[ — 0). Thus, in the results presented below we shall assume the conditions 

X > Go > Gi > G = Gi 

for the physical coupling constants. By solving Eqs. l|25H29ll . analytical solutions for the gap equations at temperature 
T = can be obtained. Because of G = G\, A v is zero for any doping x, and the SU(4) gap is 



A SU4 = ±GifVl-£ 2 - (30) 
Solution of other gaps at T = indicates that there is a critical doping point 



W^N'-f (31> 

separating qualitatively distinct solutions for the quasiparticle equations. When x > x qi there exists a pure pairing 
state [all gaps entering Eqs. l(2"5|l - (f2"§|l except the singlet pairing gap Ad vanish]. When x < x qi the gap equations 
have two solutions: the pairing solution remains but it is an excited state; another solution becomes the ground 
state, which differs from the pure pairing state by having both the AF gap A q and the pairing gap Ad nonzero, and 
becoming a pure AF state at half filling. Specifically, one finds that for x < x q 



A « = ^\l{x^-x){x q -x) (32) 



A„ 



2 

G' n 

2 
< 

''1 



\j x(x q 1 - x) (33) 
A' = -*^ Xq (l - XqX ), (34) 
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FIG. 4: A generic gap diagram for energy gaps versus doping, as predicted by the SU(4) model. The energy gaps are scaled by 
the strength of the pairing correlation G' and the doping is scaled by the maximum doping Pf. The critical point x q = 0.68 is 
chosen in this plot and x' = 17, Go = 8.2, and Gi = 9.5. 



while for < x < 1 (but this is the ground-state solution only for x > x q ), 

A, = (35) 
A = A d = ^7(1-^2) (36) 

A' = -&±x. (37) 

A generic gap diagram at T = describing features of the energy gaps as functions of the doping x is shown in Fig. 
4. The diagram is constructed using the analytical expressions ijBTjl l3"21 15"3l l3~5l l3"tJj) . It is generic because the relative 
size of the gaps can be modified by the choice of interaction strengths but their basic forms are dictated entirely by 
the algebraic structure of the model. Four energy gaps are predicted: 

1. The gap A q measuring antiferromagnetic correlations (see Eq. Q32J), which has its maximal value at x = 0, 
decreases nearly linearly to the region of the pairing gaps as doping increases, crosses the pairing gaps, and 
vanishes eventually at the critical doping x q . 

2. The spin-singlet pairing gap A^ (see Eq. (f 331) ) . which is the superconducting gap for x < x q . Starting from the 
critical doping x q , it decreases as the doping decreases. 

3. The spin-singlet pairing gap Ao = A<j (see Eq. (|36p). which is the superconducting gap for x > x q but is not 
the order parameter for the ground state in the doping range x < x q . In x < x q , it increases as the doping 
decreases. 

4. The total SU(4) correlation measured by Asu4, which has its doping dependence following Eq. (|30() in the whole 
doping range. 

It is interesting to observe that at the critical doping value x q the pairing gaps exhibit non-trivial behavior: for x > x q , 
the singlet pairing gap corresponds to a single curve, labeled by Ao. However, as the doping decreases from x q , the 
energy scale A splits into two curves (labeled Ad and A ) having very different doping dependence. 
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C. Nature of the pseudogaps in cuprates 

The qualitative features of the SU(4) gap diagram in Fig. 4 seem to agree well with many recent observations 
in cuprates^. In particular, the appearance of a quantum critical point x q , the splitting of the pairing gap in the 
underdoped regime, and the existence of two distinct pseudogaps are basic predictions of the SU(4) model that have 
some experimental support. 

The occurrence of the doping critical point and the splitting of the pairing gap below the critical point are understood 
in the SU(4) model as a direct consequence of the competing pairing and AF correlations below x q . In this picture, 
the critical doping point x q defines a natural boundary between overdoped and underdoped regimes. It corresponds 
to the disappearance of AF correlations and separates a doping regime characterized by weak superconductivity and 
reduced pair condensation energy from a doping regime characterized by strong superconductivity and maximal pair 
condensation energy. 

The origin of pseudogaps is one of the unsolved major problems in high-T c superconductor physics. The pseudogap 
is believed to have a profound effect on a wide range of physical properties of cuprates and is generally considered to 
be intimately connected with the origins of high-temperature superconductivity. As noted in Refs^, there is as yet 
little consensus about their origin and behavior but many proposed explanations fall into one of two categories: the 
preformed pair picture^ or the competing order picture^. As we now briefly discuss (more detail will be published 
elsewhereiS) , the SU(4) model predicts, and unifies at the microscopic level, two sources for pseudogap behavior, one 
of which may be interpreted as arising from preformed pairs and one of which may be interpreted as arising from 
competing order. 

The preformed pair picture^ suggests that the pseudogap is the energy that is needed to form pairs as a first step 
in condensing them into a state with long-range order. A pseudogap arising in this way would be expected to decrease 
with increased doping and should merge with the pairing gap in the overdoped region. The SU(4) gap Agu(4) has 
exactly these properties. As can be seen from Eq. (|3(J[1 . Agu(4) depends only on G\ and is independent of x q . It is 
maximal at half-filling and falls to zero at the largest doping. This gap is related to the SU(4) invariant £ SU 4 (see Eq. 
(1191) ). which is state independent. Therefore, the SU(4) gap is expected to have no temperature dependence, as long 
as the SU(4) symmetry is preserved. 

On the other hand, the competing order picture^ suggests that the pseudogap is the energy scale for some form of 
order competing with superconductivity. A pseudogap originating in competing order should not generally merge with 
the pairing gap and may fall to zero at a critical doping point where the competing order is completely suppressed 
relative to the superconducting order. The SU(4) antiferromagnetic gap A g has these properties and gives rise to a 
pseudogap that may be interpreted in terms of competing AF and SC order. Specifically, the scale A 9 is the energy 
per electron required to break the AF correlation and serves as an order parameter of the AF phase. 

Thus, the SU(4) model predicts the existence of two distinct pseudogaps. The seemingly rather contradictory 
mechanisms for the origin of the pseudogaps as suggested in Refsi2^2& may find a consistent explanation in the SU(4) 
model. 



VI. CONCLUSION 

The SU(4) symmetry represents a fully microscopic fcrmion system in which SC and AF modes enter on an equal 
footing. At this "unification" level, there is in a sense no distinction among these degrees of freedom, just as in the 
Standard Electroweak Theory of elementary particle physics the electromagnetic and weak interactions are unified 
above the intermediate vector boson mass scale. We may view the system as having condensed into SU(4) pairs, 
which fixes the length of the state vectors (through the SU(4) Casimir expectation value) but not their direction in 
the state-vector space. Physically, this means that the system is paired, with the pair structure exhibiting SU(4) 
symmetry, but is neither SC nor AF because fluctuations on a scale set by the temperature prevent selection of 
the SC or AF directions. Stated in another way, the SU(4) pairs are of collective strength, but are not condensed 
into a state with long-range order. Stated in yet another way, neither the AF nor SC order parameters have finite 
expectation values in this regime, but a sum of their squares (the SU(4) Casimir) does. This constraint implies an 
intimate connection between superconductivity and antiferromagnetism in the SU(4) model. They are, in a sense, 
different projections of the same fundamental vector in an abstract algebraic space. 

Compared to the Hubbard or t-J models, the dynamical symmetry approach applied here represents a different 
way to simplify a strongly-correlated many-body problem. In the Hubbard or t-J models, approximations are made 
to simplify the Hamiltonian but no specific truncation is assumed for the configuration space, although in practical 
calculations a truncation is typically necessary. In contrast, we make no approximation to the Hamiltonian. The only 
approximation is the (severe) space truncation. The symmetry dictated Hamiltonian includes all possible interactions 
in the truncated space. In principle, if all the degrees of freedom of the system are included, this approach constitutes 
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an exact microscopic theory. The validity of this approach depends entirely on the validity of the choice of truncated 
space, which may be tested by comparison with the data. 

The advantage of the dynamical symmetry approach is its cleanness and simplicity. It is clean because the only 
approximation is the selection of the truncated model space. Thus, a failure of the method is a strong signal that one 
has chosen a poor model space. It is simple because the method supplies analytical solutions for various dynamical 
symmetry limits as a starting point. These symmetry- limit solutions provide an immediate handle on the physics and 
permit an initial judgement of the model's validity without large-scale numerical calculations. Beyond the symmetry 
limits, numerical calculations may be necessary. However, because of the low dimensionality of the models spaces and 
the power of group theory, such numerical calculations are much simpler than those in, say, a Hubbard or t-J model. 

Simple symmetries as a predictor of dynamics has found powerful application in fields such as particle physics or 
nuclear physics. Although a real physical system may be very complex, these applications are based on the obvious 
point that nature has managed to construct a stable ground state having well-defined collective properties that change 
systematically from compound to compound. Extensive experience in many fields of physics suggests that this is the 
signal that the phenomenon in question is described by a small effective subspace with renormalized interactions (that 
may differ substantially in form and strength from those of the full space) and governed by a symmetry structure of 
manageable dimensionality. The results described here suggest that the observational complexities of high temperature 
superconductors may also be described by a simple theory based on dynamical symmetries having these properties. 
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The SU(4) model is particle-hole symmetric. For cuprate superconductors one is generally interested in hole-doped com- 
pounds but more general applications of SU(4) may be interested in either electrons or electron holes. Hereafter, unless 
specified explicitly, we use "electrons" to reference either electrons or holes. 

Groups generally may have more than one Casimir invariant. We shall use the term "Casimir" to refer loosely to the lowest- 
order such invariants (which are generally quadratic in the group generators). In the context of the present discussion, 
quadratic Casimirs are associated with 2-body interactions at the microscopic level. Higher-order Casimirs are then generally 
associated with 3-body and higher interactions. The restriction of our Hamiltonians to polynomials of order 2 in the Casimirs 
is then a physical restriction to consideration of only 1-body and 2-body interactions. 



